Lifshitz Interaction between Dielectric Bodies of Arbitrary Geometry 
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A formulation is developed for the calculation of the electromagnetic-fluctuation forces for dielec- 
tric objects of arbitrary geometry at small separations, as a perturbative expansion in the dielectric 
contrast. The resulting Lifshitz energy automatically takes on the form of a series expansion of the 
different many-body contributions. The formulation has the advantage that the divergent contribu- 
tions can be readily determined and subtracted off, and thus makes a convenient scheme for realistic 
numerical calculations, which could be useful in designing nano-scale mechanical devices. 
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Since the pioneering work of van der Waals in 1873 — 
which revealed that condensation of gases is not possible 
unless an attractive interaction is at work in matter — and 
the explanation of this interaction by London in 1930 in 
terms of the correlation between quantum fluctuations 
of atoms and molecules [l[, we have been witnessing a 
surge of interest in exploring the physical implications 
of such fluctuation-induced interactions @, Q • The case 
of two parallel plates made of perfect conductors tackled 
by Casimir 0] and the subsequent generalization to the 
case of dielectric materials by Lifshitz 5i] made important 
contributions to our understanding of the macroscopic 
manifestations of these interactions. Relevant experi- 
mental studies, which were ongoing alongside with the 
theoretical developments [lj, culminated recently with 
high precision quantitative verifications of the Casimir 
force BR SI- 

The advent of nanotechnology in recent years has 
added to the interest in electromagnetic-fluctuation in- 
teractions as they become dominant at nanoscale [Io| . 
and quantitative knowledge of them appears to be nec- 
essary in designing nano-machines to avoid unwanted 



effects such as stiction [ill, ll2(. Moreover, one could 
even make use of these interactions, e.g. the normal 
Casimir force between a flat plate and a sphere [1, [l3j 
or the lateral Casimir force between two corrugated sur- 
faces 14, Il5j |. in designing novel actuation schemes in 



microelectromechanical systems (MEMS). It thus seems 
desirable to be able to calculate the electromagnetic- 
fluctuation forces for a given assortment of dielectric and 
metallic objects in a certain configuration. 

This, however, appears to be a nontrivial task due a 
number of subtleties involved. Casimir forces are known 
to depend in a nontrivial way on the geometry of the ob- 
jects [16], and the theoretical schemes that have so far 
been developed to study this effect are only applicable to 
perfect conductors 0, Hi, El H EI H3| ■ On the other 
hand, Lifshitz has pointed out that the electromagnetic- 
fluctuation force between two boundaries is dominated 
by the dielectric properties of the media in the frequency 
that is set by the separation between them This 
means that when these forces are most relevant, i.e. at 
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FIG. 1: Schematics of (a) the dielectric function profile, 
and the many-body contributions to the electromagnetic- 
fluctuation forces, and (b) the two semi-infinite dielectric bod- 
ies with irregular surfaces. 



length scales lower than 100 nm that is of the order the 
plasma wavelengths of most good metals, the perfect con- 
ductor assumption in the calculation of the Casimir force 
breaks down. Finally, the dependence of the divergent 
contributions to the Casimir energy — that should be re- 
moved in a carefully regularized formulation — on the ge- 
ometry of the boundaries is not well characterized, and 
this makes it difficult to develop systematic numerical 
schemes of calculations. 

The fact that the Lifshitz interactions at small separa- 
tions effectively involve dielectric constants at relatively 
high-frequencies suggests that a strategy based on expan- 
sion in dielectric contrast could act as a useful comple- 
mentary approach to the existing formulations. Here, we 
have developed a path integral formulation for calculat- 
ing the Lifshitz energy for dielectric bodies of arbitrary 
geometry, as an expansion in powers of the variations in 
the dielectric constant in space (see Fig. QJ,). The result 
effectively takes on the form of an expansion in many- 
body interactions with a fundamental tensorial kernel 
that reflects the nature of the electromagnetic fluctua- 
tions, with explicit expressions for all of the terms in the 
series. It has the advantage that the divergent contribu- 
tions in the series are manifest and can be subtracted off 
systematically, which makes it very suitable for numerical 
calculations. 

We start with the action integral for the electromag- 
netic field in a matter with dielectric constant e, i.e. 
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J dtd 3 x [eE 2 — B 2 1 , where standard expressions 
for the electric field E = —-dtA — V<p and the magnetic 
field B = VxA are assumed in terms of the potentials 
(j) and A. To perform the quantization, we need to ex- 
press the action in terms of the potentials and break the 
gauge invariance by choosing a particular gauge. Using 
the temporal gauge, cj> = 0, we can write the action inte- 
gral as S m = ± J <g J <p T Ai (-u>, r)[e(a,, r)( W /c) 2 % + 
dkdkSij — didj]Aj (co,r), where summation over repeated 
indices is assumed. In this expression, e(w,r) is the 
frequency dependent dielectric function of the medium, 
which could represent any spatial arrangement of dielec- 
tric objects of arbitrary shapes, as depicted in Fig. 

The quantization can now be performed by using the 
path integral method, which is facilitated if a Wick rota- 
tion is performed in the frequency domain. This renders 
the action integral "Euclidean" and we can find the par- 
tition function Z = J VA % exp (-S® uclidcan /h), and the 
Lifshitz energy E = — y InZ, where T is a long observa- 
tion time. The calculation yields 
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where K-i 

r'). Writing e(i(,r) = 1 + 5e(i(,r), we can decom- 
pose the kernel Kij into a diagonal part that corre- 
sponds to the empty space and a perturbation that en- 
tails the dielectric inhomogeneity profile. In Fourier 
space, this reads /Cjj(C;q,q') = /Co,y(C, q)(27r) 3 5 3 (q + 
q') + SJCijiC^q 1 ), where /C ,.y(C,q) = (C/ c ) 2 ^ + 



q t qj and SJdj {(; q, q') = ((/c) 2 S lj Se(iC, q 



q'). We can now recast the expression for the 
Lifshitz energy into a perturbative series by using 
trln[/C] = trln[£ ] + £"=i tr[(/C^/C)"], where 

(C/cfdij + q l q j 



^n«(C ! q) = 



(C/c) 2 [(C/ C ) 2 



Using these definitions, 



one can write the explicit form for the trace as 

[(C/crS lll2 + q^q^} ■ ■ ■ [(C/cf6 inil + $>qW] 



[(C/C) 2 + <J« 2 ] • • • [(C/C) 2 + 9 (») 2 



x 5e(iC,-q (1) +q (2) )--- 6e(i(, 



, (2) 



which involves the geometric information about the ar- 
rangement of the dielectric objects through the Fourier 
transform of the dielectric function profile. Transform- 
ing back the expression in Eq. @ into real space, we 
find the following series for the Lifshitz energy of any 
heterogeneous dielectric medium 
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FIG. 2: Comparison of the expansion in Eq. (|3} up to the sec- 
ond, fourth, and sixth order, as well as the Clausius-Mossotti 
approximation, with the exact Lifshitz result for a pair of par- 
allel semi-infinite dielectrics, as a function of their separation 
H in units of the plasma wavelength A p . The dielectric con- 
stant of both objects is chosen as e(i£) = 1 + ui^,/(uiq + £ 2 ) 
with ujq = 2 Wp. 



where the Green's function associated with the electro- 
magnetic fluctuations is defined as 
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This result has a number of interesting characteristics. 
First, it appears that an expansion in powers of Se auto- 
matically turns into a summation of integrated contribu- 
tions of n-body interactions. Moreover, all of the n-body 
interaction terms have simple expressions in terms of a 
single fundamental kernel 0«(r) that mediates the two- 
body part of the interaction. This kernel is tensorial, and 
has the structure of the electric field of a radiating dipole 
in imaginary frequency 211 ] . It has, in fact, been intro- 
duced some time ago in connection with van der Waals 
interactions [12] ■ Finally, the result has a closed form 
expression for the Lifshitz energy for any geometrical ar- 
rangement of dielectric bodies in terms of quadratures. 

To examine the convergence property of the above 
series, we can take a specific example for which the 
exact result is known, and compare it with what 
we find from the first few terms in Eq. ([3]). Wc 
consider the case of two identical semi-infinite dielec- 
tric objects that are placed parallel to each other 
at a separation H, for which the exact expression 
for energy per unit area is known to be ^Lif = 
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(3) To calculate the energies, we have assumed a simple form 



| , where s = \J e — 1 
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for the dielectric constant, 



where uj p represents the plasma frequency, from which 
the plasma wavelength A p = 2ttc/lu p can be extracted. 
In Fig. [21 the ratio between the energy as calculated 
from Eq. §3§ up to the second, fourth, and sixth or- 
der and the exact Lifshitz result is shown as a function 
of the separation in units of the plasma wavelength for 
luq = 2 u> 2 and 7 = 0, which corresponds to a dielectric 
with e(0) = 1.5. The results clearly show a crossover 
between two asymptotic regimes near H « A p /(27r) [23j |. 
The series appears to be rapidly convergent so long as 
5e(0) < 1, and the convergence is considerably more ef- 
ficient for H < X p /(2n) [as compared to H > A p /(27r)], 
especially for <5e(0) < 1. The damping constant 7 ap- 
pears to have very little effect on the Lifshitz energy for 
dielectrics (i.e. when w 7^ 0). The results presented 
in Fig. [21 for example, do not change appreciably for 
nonvanishing damping constants of up to 7 = 0.1 uj p . 

It is quite well known that any calculation of the 
Casimir or Lifshitz energy encounters a variety of di- 
vergent contributions that are very difficult to charac- 
terize. For example, it is not clear how these divergent 
terms depend on the geometry of the system, so that they 
could be identified in a simple geometry and done away 
with in a systematic way for slightly deformed bound- 
aries, in the calculations that involve perturbation in the 
geometry of the objects 14, l]| . While the present for- 
mulation also suffers from this deficiency, in the sense 
that the expression in Eq. ^ involves divergent con- 
tributions, the fact that the expressions for the vari- 
ous terms in the series are known a priori allows for 
systematic identification of these contributions and thus 
their systematic cancelation. For example, one can show 
that putting 1 instead of the fraction in the second line 
of Eq. @ yields a singular contribution of the form 

Ssi„ g . = fiff/A In[e(iC,r)]/^r, which could 
be subtracted off systematically (i.e. at every order n 
in the series expansion) for any geometry and dielectric 
configuration. 

Let us now focus our attention on the specific arrange- 
ment shown in Fig. [TJ), where two semi-infinite dielec- 
tric bodies with irregularly shaped boundaries are placed 
nearly parallel to each other at a mean separation H . We 
can write down the dielectric function profile as 

' e 2 «), H + /l 2 (x) < z < +00, 

e(iC,r) = <l, fti(x) < z <H + h 2 (x), (5) 

k ei«), -00 < z < fci(x), 

and the corresponding Fourier transform as 5e(iC, q) = 
£ / d 2 x e lcu ' x [<5e 2 e lq ^ H+h2 ^ - 8ei e iq * hl W]. When 
the separation of the surfaces is smaller than the plasma 
wavelengths of the two dielectric media, and Se± and Se 2 



are small compared to unity, the leading contribution in 
Eq. ([3]) comes from the second order term. Putting in 
the dielectric function profile, we find 
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with the incomplete gamma function defined as T(a, z) 



a-l 



Note that at this order, the Lifshitz 



energy is pairwise additive. In the above result, we have 
kept the frequency dependence of the dielectric functions 
as well as the geometry of the boundaries arbitrary for 
generality of the presentation [24j . The expression in Eq. 
([6|) can be considerably simplified for /ii(x) = 0: 



£ 2 | Mx)=0 = J d 2 x£ 2 (H + h 2 ( X j) 



(8) 



in terms of the original Lifshitz result for the energy per 
unit area of flat boundaries [B| 



64tt 2 c 2 



(9) 

with C{u) = A(u 2 - l)£a(2u) + + 2u + u 2 - 2w 3 ), 

and the exponential integral function defined as E n (z) — 
J. dt e~ zt /t n . This simplification is a general feature of 
pairwise additive interactions, as shown in Ref. (l8j . 

The pairwise summation approximation is widely used 
in the literature not only for dielectric materials at close 
separations but also for perfect conductors at arbitrary 
separations. We can use the present formulation to shed 
some light on the nature of this approximation, and make 
assessment on what is not included. Going back to Eq. 
(|31), we note that the Green's function S^-(r) is decom- 
posed into a long-ranged kernel, and a "contact" part in 
the form of |<%<$ (r). A systematic pairwise summation 



approximation then amounts to keeping two of the ^-'s 
in its full form and approximating the rest of them by 
their contact-contribution, in each term of the series in 
Eq. ([3]). Considering the different ways of doing this and 
keeping track of the corresponding combinatorial prefac- 
tors, one can then sum up the entire series and find a 
closed form expression for the Lifshitz energy between 
objects of arbitrary shapes in the pairwise summation 
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approximation. The result will be identical to Eq. ([6]) ex- 
cept for the replacement <5ej — > 1+ ^ £ e '.y 3 - ) = 3 ( j-^ ^ ■ The 
combination reminds us of the Clausius-Mossotti equa- 
tion for molecular polarizability, which is well known to 
be valid only for dilute materials such as gases [21] . This 
approximation allows us to calculate the energy for per- 
fect conductors, and we can use the example of two flat 
and parallel perfect conductors, for which the exact result 
is known to be f^xact = ~ TMTP"' to assess its validity. 
Letting ei,e 2 -> oo, we find £p WS = - ga ; which 
compares to the exact result as fe WS = l^r — 0.797. 
While this ratio varies with geometry, it gives us an es- 
timate of typical errors that are involved in the calcula- 
tions based on pairwise summation approximation. We 
also note that any attempt in going beyond this approxi- 
mation should include the tensorial structure involved in 
Eq. |(3J|, and for example, an ad hoc augmentation by in- 
troduction of scalar three-body interactions etc. will not 
be justified in light of our present scheme. The Clausius- 
Mossotti approximation appears to give comparatively 
better results for the case of dielectrics, as can be seen 
from the example in Fig. O 

The formulation presented here could also be ap- 
plied to the case of magnetic materials j25[. For a 
medium that is described by the dielectric function pro- 
file e(w,r) and the magnetic permeability profile /i(w,r), 
we should replace the kernel in Eq. ([T]) by /C«(£; r, r') = 
[(C/c) 2 e(*C r)% + di^dj ~ d k7 ^d k 6 ZJ ]S 3 (r - r'), 
and the rest of the procedure follows closely. Generaliza- 
tion to the case of finite temperatures by discretizing the 
frequency is also straightforward. 

In conclusion, we have presented a path integral for- 
mulation for the calculation of the Lifshitz energy for 
dielectric materials of arbitrary shape, as a series ex- 
pansion in the dielectric contrast. The expansion con- 
verges very rapidly for dielectric objects that are at sep- 
arations considerably smaller than their corresponding 
plasma wavelengths, and is expected to work perfectly 
for surfaces at nanometric separations. The results pre- 
sented here could be applicable for the calculation of 
electromagnetic-fluctuation forces that are involved in 
nano-mechanical devices. 

It is a great pleasure to acknowledge fruitful discussions 
with M.A. Charsooghi, T. Emig, A. Hanke, M. Kardar, 
R. Matloob, M.F. Miri, and F. Mohammad-Ranee. 
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